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Performance analysis of spatial smoothing 
schemes in the context of large arrays. 

Gia-Thuy Pham, Philippe Louhaton, Fellow, IEEE and Pascal Vallet, Member, IEEE, 


Abstract 

This paper adresses the statistical behaviour of spatial smoothing subspace DoA estimation schemes 
using a sensor array in the case where the number of observations N is significantly smaller than the 
number of sensors M, and that the smoothing parameter L is such that M and NL are of the same order 
of magnitude. This context is modelled by an asymptotic regime in which NL and M both converge 
towards cxd at the same rate. As in recent works devoted to the study of (unsmoothed) subspace methods 
in the case where M and N are of the same order of magnitude, it is shown that it is still possible 
to derive improved DoA estimators termed as Generalized-MUSIC with spatial smoothing (G-MUSIC 
SS). The key ingredient of this work is a technical result showing that the largest singular values and 
corresponding singular vectors of low rank deterministic perturbation of certain Gaussian block-Hankel 
large random matrices behave as if the entries of the latter random matrices were independent identically 
distributed. This allows to conclude that when the number of sources and their DoA do not scale with 
M, N, L, a situation modelling widely spaced DoA scenarios, then both traditional and Generalized 
spatial smoothing subspace methods provide consistent DoA estimators whose convergence speed is 
faster than The case of DoA that are spaced of the order of a beamwidth, which models closely 
spaced sources, is also considered. It is shown that the convergence speed of G-MUSIC SS estimates is 
unchanged, but that it is no longer the case for MUSIC SS ones. 

I. Introduction 

The statistical analysis of subspace DoA estimation methods using an array of sensors is a topic that 
has received a lot of attention since the seventies. Most of the works were devoted to the case where the 
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number of available samples N of the observed signal is much larger than the number of sensors M of the 
array (see e.g. ifldll and the references therein). More recently, the case where M and N are large and of 
the same order of magnitude was addressed for the first time in lITOl using large random matrix theory. lITOl 
was followed by various works such as Q, ifTTl . lUl, IH. The number of observations may also be much 
smaller than the number of sensors. In this context, it is well established that spatial smoothing schemes, 
originally developped to address coherent sources (||4l, |[T3l . ifTTI l. can be used to artificially increase the 
number of snapshots (see e.g. |[T4l and the references therein, see also the recent related contributions ifTSl . 
|[T^ devoted to the case where N = 1). Spatial smoothing consists in considering L < M overlapping 
arrays with M — L + 1 sensors, and allows to generate artificially NL snapshots observed on a virtual 
array of M — L + 1 sensors. The corresponding (M — L + 1) x A^L matrix, denoted collecting the 

observations is the sum of a low rank component generated by (M —L + l)-dimensional steering vectors 
with a noise matrix having a block-Hankel structure. Subspace methods can still be developed, but the 
statistical analysis of the corresponding DoA estimators was addressed in the standard regime where 
M — L + 1 remains fixed while NL converges towards oo. This context is not the most relevant when 
M is large because L must be chosen in such a way that the number of virtual sensors M — L + 1 be 
small enough w.r.t. NL, thus limiting the statistical performance of the estimates. In this paper, we study 
the statistical performance of spatial smoothing subspace DoA estimators in asymptotic regimes where 
M — L + 1 and NL both converge towards oo at the same rate, where ^ 0 in order to not affect the 

aperture of the virtual array, and where the number of sources K does not scale with M, N, L. For this, 
it is necessary to evaluate the behaviour of the K largest eigenvalues and corresponding eigenvectors of 
the empirical covariance matrix —■ To address this issue, we prove that the above eigenvalues 

and eigenvectors have the same asymptotic behaviour as if the noise contribution to matrix 
a block-Hankel random matrix, was a Gaussian random matrix with independent identically distributed. 
To establish this result, we rely on the recent result HI addressing the behaviour of the singular values of 
large block-Hankel random matrices built from i.i.d. Gaussian sequences. HI implies that the empirical 
eigenvalue distribution of matrix ^ — converges towards the Marcenko-Pastur distribution, and that 
its eigenvalues are almost surely located in the neigborhood of the support of the above distribution. This 
allows to generalize the results of HI to our random matrix model, and to characterize the behaviour of 
the largest eigenvalues and eigenvectors of —. We deduce from this improved subspace estimators, 

called DoA G-MUSIC SS (spatial smoothing) estimators, which are similar to those of ifTTIl and HI. We 
deduce from the results of mi that when the DoAs do not scale with M, N, L, i.e. if the DoAs are 
widely spaced compared to aperture array, then both G-MUSIC SS and traditional MUSIC SS estimators 
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are consistent and converge at a rate faster than . Moreover, when the DoAs are spaced of the order of 
the behaviour of G-MUSIC SS estimates remains unchanged, hut the convergence rate of traditional 
suhspace estimates is lower. 

This paper is organized as follows. In section ini we precise the signal models, the underlying as¬ 
sumptions, and formulate our main results. In section JIIJ we prove that the largest singular values and 
corresponding singular vectors of low rank deterministic perturbation of certain Gaussian block-Hankel 
large random matrices behave as if the entries of the latter random matrices were independent identically 
distributed. In section HVl we apply the results of section ITm to matrix and follow Q in order 

to propose a G-MUSIC algorithm to the spatial smoothing context of this paper. The consistency and 
the convergence speed of the G-MUSIC SS estimates and of the traditional MUSIC SS estimates are 
then deduced from the results of ifTSl . Finally, section |V] present numerical experiments sustaining our 
theoretical results. 

Notations : For a complex matrix A, we denote by A^, A* its transpose and its conjugate transpose, 
and by Tr (A) and ||A|| its trace and spectral norm. The identity matrix will be I and will refer to 
a vector having all its components equal to 0 except the n-th equals to 1. For a sequence of random 
variables (A„)„gN and a random variable X, we write 

n—>-cxD 

when Xn converges almost surely towards X. Finally, = Op (1) will stand for the convergence of A„ 
to 0 in probability, and X^ = Op{l) will stand for tightness (boundedness in probability). 

II. Problem formulation and main results. 

A. Problem formulation. 

We assume that K narrow-band and far-held source signals are impinging on a uniform linear array 
of M sensors, with K < M. In this context, the M-dimensional received signal (yn)n>i can be written 
as 

Yn — T Wjj, 

where 

• Am = [slm{0i), ..., slm{0k)] is the MxK matrix of M-dimensionals steering vectors a.M{0i), ■ ■ ■, slm{0k), 

with 6i,... ,9 k the source signals DoA, and aLM{9) = ; 

• s„ S contains the source signals received at time n, considered as unknown deterministic ; 
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• (vn)n>i is a temporally and spatially white complex Gaussian noise with spatial covariance E[v„v* ] = 

The received signal is observed between time 1 and time N, and we collect the available observations 
in the M x N matrix Yn defined 


Yat = [yi,... ,y7v] = AmSat + Vtv, 


( 1 ) 


with Sat = [si,..., sat] and Vat = [vi,..., vat]. We assume that Rank(SAr) = K for each M, N greater 
than K. The DoA estimation problem consists in estimating the K DoA 9i,... ,9 k from the matrix of 
samples Yat. 

When the number of observations N is much less than the number of sensors M, the standard subspace 
method fails. In this case, it is standard to use spatial smoothing schemes in order to artificially increase 
fhe number of observafions. In particular, if is well esfablished fhaf spatial smoofhing schemes allow fo 
use subspace mefhods even in fhe single snapshof case, i.e. when N = 1 (see e.g. |[T4l and fhe references 
fherein). If L < M, spatial smoofhing consisfs in considering L overlapping subarrays of dimension 
M — L + 1. Af each time n, L snapshofs of dimension M — L + 1 are fhus available, and fhe scheme 
provides NL observations of dimension M — L + 1. In order fo be more specific, we infroduce fhe 
following nofafions. If L is an infeger less fhan M, we denote by fhe (M — L + 1) x L Hankel 
matrix defined by 

yi,n y2,n . y L,n 



y2,n y3,n 


y L-\-l,n 


( 2 ) 


\ yM-L+l,n yM-L+2,n . yM,n / 

Column I of mafrix corresponds fo fhe observation on subarray I af time n. Collecfing all fhe 
observations on fhe various subarrays allows fo obfain NL snapshofs, fhus increasing artificially fhe 
number of observafions. We define Y^^ as fhe (M — L + 1) x NL block-Hankel mafrix given by 


^(L) _ 
^ N — 



( 3 ) 


In order fo express Y^j^\ we consider fhe (M — L + 1) x L Hankel mafrix A^^\9) defined from vector 
^M{d) in the same way fhan y^\ We remark fhaf A^^\9) is rank 1, and can be written as 

A^^\9) = ^L{M -L + 1)/M aM-L+i(0) (4) 
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We consider the (M — L + 1) x KL matrix 

aW = .. .,A^^\eK)) (5) 

which, of course, is a rank K matrix whose range coincides with the suhspace generated hy the (M — 
L + 1)-dimensional vectors aLM-L+i{di), ■ ■ ■, can he written as 

) = A(^) (Sat ® II) + VSJ) (6) 


where matrix is the block Hankel matrix corresponding to the additive noise. As matrix Sat^Il is 
full rank, the extended ohervation matrix y]^^ appears as a noisy version of a low rank component whose 
range is the A-dimensional suhspace generated hy vectors aLM-L+i{Gi), ■ ■ ■, aAf_x,+i(0A')- Moreover, 
it is easy to check that 

E 


'vWyW*' 

NL 


= aA 


M-L+l 


Therefore, it is potentially possible to estimate the DoAs {Ok)k=i,...,K using a subspace approach based 
on the eigenvalues / eigenvectors decomposition of matrix Yj^^Y^^*/AL. The asymptotic behaviour of 
spatial smoothing subspace methods is standard in the regimes where M — L + 1 remains fixed while NL 
converges towards oo. This is due to the law of large numbers which implies that the empirical covari¬ 
ance matrix /NL has the same asymptotic behaviour than A^^^ (SatS^ (8* Il/NL) 

a‘^lM-L+ 1 , In this context, the orthogonal projection matrix onto the eigenspace associated to the 
M—L+l—K smallest eigenvalues of /NL is a consistent estimate of the orthogonal projection 

matrix 11^^^ on the noise subspace, i.e. the orthogonal complement of sp{aAf-L-i-i(6*i)) ■ ■ ■ j ^m-l+i/^k)}- 
In other words, it holds that 




0 a.s. 


(7) 


The traditional pseudo-spectrum estimate {6) defined by 


flt\e) = SiM-L+lieyiVN’^M-L+m 


(L) 


fhus verifies 


sup 

0e[—7r,7r] 




N^oo 


-> 0 . 


( 8 ) 


where ri{6) = a.M-L+iiG)*Ll^^^a.M-L+i{G) is fhe MUSIC pseudo-specfrum. Moreover, fhe K MUSIC 
fradifional DoA esfimafes, defined formally, for A; = 1,..., A, by 


= argmin57jv^(6l), 
O^Xk 




(9) 
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where is a compact interval containing 9^ and such that n X; = 0 for A; / I, are consistent, i.e. 


'k,N 


N^oo 


> Ok- 


( 10 ) 


However, the regime where M — L + 1 remains fixed while NL converges towards oo is not very 
interesting in practice because the size M — X + 1 of the suharrays may he much smaller that the number 
of antennas M, thus reducing the resolution of the method. We therefore study spatial smoothing schemes 
in regimes where the dimensions M — X + 1 and NL of matrix are of the same order of magnitude 
and where ^ 0 in order to keep unchanged the aperture of the array. More precisely, we assume that 

integers N and X depend on M and that 

, 1 M - X + 1 

M^+oo,iV = 0(M^), - </3< 1, cjv =--(11) 

In regime (fTTl) . N thus converges towards oo but at a rate that may be much lower than M thus 
modelling contexts in which N is much smaller than M. As N ^ +oo, it also holds that ^ c*. 
Therefore, it is clear that X = 0{M°) where a = 1 — (3 verifies wifh 0 < a < 2/3. X may fhus converge 
towards oo (even faster fhan if /3 < 1/2) buf in such a way fhat ^ 0. As in regime (fTTl) N depends 

on M, it could be appropriate to index the various matrices and DoA estimators by integer M rather 
than by integer N as in definitions (l5]l and ®. However, we prefer to use the index N in the following 
in order to keep the notations unchanged. We also denote projection matrix and pseudo-spectrum 
r/(0) by and rji^{6) because they depend on M. Moreover, in the following, the notation N —)■ +oo 
should be understood as regime (fTTl) for some f3 G (1/3,1]. 


B. Main results. 

In regime (fTTl) . (fT]) is no more valid. Hence, (ITOl) is questionable. In this paper, we show that it is 
possible to generalize the G-MUSIC estimators introduced in Q in the case where X = 1 to the context 
of spatial smoothing schemes in regime (fTTl) and establish the following results. Under the separation 
condition that the K non zero eigenvalues of matrix are above the threshold 

for each N large enough, we deduce from ifTSl that: 

• the spatial smoothing traditional MUSIC estimates {0^'^jKf)k=i,...,K and the G-MUSIC SS estimates, 
denoted {6k,N)k=i,...,K are consistent and verify 

M(0g}v-0fc) ^ Oa.s., (12) 

M{9k,N — Ok) —)• 0 a.s. (13) 
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(fT^ and (fT3]) hold when the DoA {6k)k=i,...,K do not scale with M,N. In pratice, this assumption 
corresponds to practical situations where the DoA are widely spaced because when the DoA {9k)k=i,...,K 
are fixed, the ratio 

\0k-0i\ 

( 27 r)/M 

converges towards oo. We deduce from ifTSl that: 

• If A' = 2 and that the 2 Do As scale with M, N is such a way that 62 ,n — 0i,N = then the 

G-MUSIC SS estimates still verify (fTSl) while the traditional MUSIC SS estimates no longer verify 

CHI 

As in the case L = 1, the above mentioned separation condition ensures that the K largest eigenvalues 

of the empirical covariance matrix correspond to the K sources, and the signal and 

noise subspaces can be separated. In order to obtain some insights on this condition, and on the potential 

benefit of the spatial smoothing, we study the separation condition when M and N converge towards 00 

at the same rate, i.e. when ^ d^, or equivalently that j3 = 1 and that L does not scale with N. In this 

s s* 

case, it is clear that c* coincides with c* = d^/L. Under the assumption that converges towards a 

diagonal matrix D when N increases, then we establish that the separation condition holds if 

\k ( 14 ) 

for each (M, N) large enough. If L = 1, the separation condition introduced in the context of (un¬ 
smoothed) G-MUSIC algorithms (Q) is of course recovered, i.e. 

(A^AmD) > 


If M is large and that L « M, matrix A^_j^_|_^ Am-l-i-i is close from A\^Am and the separation 
condition is nearly equivalent to 


Ax (A^AmD) > 


O" 


'y/d^ 


Vl 


Therefore, it is seen that the use of the spatial smoothing scheme allows to reduce the threshold a'^y/d^ 
corresponding to G-MUSIC method without spatial smoothing by the factor y/L. Therefore, if M and 
N are the same order of magnitude, our asymptotic analysis allows to predict an improvement of the 
performance of the G-MUSIC SS methods when L increases provided L « M. If L becomes too large, 
the above rough analysis is no more justified and fhe impact of the diminution of the number of antennas 
becomes dominant, and the performance tends to decrease. 
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III. Asymptotic behaviour of the largest singular values and corresponding 


SINGULAR VECTORS OF FINITE RANK PERTURBATIONS OF CERTAIN LARGE RANDOM 

BLOCK-HANKEL MATRICES. 


In this section, N, M, L still satisfy (fTTl) while AT is a fixed integer that does not scale with N. We 
consider the (M + L — 1) x NL hlock-Hankel random matrix defined previously, and introduce 
matrix Z^r defined 


Zn = 


=V 


(L) 


in order fo simplify fhe nofations. The entries of Z^r have of course variance a'^/NL. In the following, 
Btv represents a deterministic (M + L — 1) x NL matrix verifying 


sup ||Bjv|| < +CX), Rank(BAr) = K, (15) 

N 

for each N large enough. We denote hy > \ 2 ,n ■ ■ ■ > Xk,n the non zero eigenvalues of matrix 
BatB^ arranged in decreasing order, and hy {uk,N)k=i,...,K and {'vik,N)k=i,...,K the associated left and 

right singular vectors of B^r. The singular value decomposition of Bat is thus given hy 

K 

k=l 

Moreover, we assume that: 


Assumption 1. The K non zero eigenvalues {Xk,N)k=i,...,K of matrix BtvB^ converge towards Ai > 
A 2 > . ■ ■ > \k when N + 00 . 

Here, for ease of exposition, we assume that the eigenvalues {\k,N)k=i,...,K have multiplicity 1 and 
that Afc 7 ^ A; for k 1. However, the forthcoming results can he easily adapted if some A^ coincide. 

We define mafrix Xat as 

Xat = Bat + Zat (16) 

Xat can thus he interpreted as a rank K perturbation of the hlock-Hankel matrix Z^. The purpose of this 
section is to study the behaviour of the K largest eigenvalues {Xk,N)k=i,...,K of matrix XatX^ as well 
as of their corresponding eigenvectors {iik,N)k=i,...,K- It turns out that {Xk,N)k=i,...,K and {iik,N)k=i,...,K 
behave as if the entries of matrix Zn where i.i.d. To see this, we have first to precise the behaviour of 
the eigenvalues of matrix ZatZ^ in the asymptotic regime (fTTT) . 
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A. Behaviour of the eigenvalues of matrix ZjyZ^. 

We first recall the definition of the Marcenko-Pastur distrihution of parameters cr^ and c (see e.g. 
ID). /iCT 2 c is the prohahility distrihution defined hy 


dfla^^cix) = (5o[l - C ]+ + 


_ll ^ s/{x-X ){x+ -x) 


^x-,x+]{x) dx 


2cr^C7rx 

with x~ = (T^(l — \fc)‘^ and x'^ = (7^(1 + v^)^. Its Stieltjes transform ma^^c{.z) defined hy 


is known to satisfy the fundamental equation 

= 

or equivalently, 

rria^A^) = 


X — z 


—z + cr'^T-, —2-^^- 

l+cr2cm^2_A^) 


1 


Alcr2,c(^) = 


-z{l + ,.{z)) 

1 


(17) 

(18) 
(19) 


—z{l + a’^cnia'^^ciz)) 

where mo -2 ^(z) is known to coincide with Stieltjes transform of the Marcenko-Pastur distrihution tt(T^c,c-^ = 
CHa^,c + (1 - c)5o- 

In order to simplify the notations, we denote hy m^{z) and m^{z) the Stieltjes transforms of Marcenko- 
Pastur distrihutions /Xcr^.c. and m^{z) and m^{z) verify Equations (fT^ and (fT^ for c = c*. 

We also denote hy xf and xf the terms xf = fT^(l — and xf = a‘^(l + We recall that 

function tu* (z) defined hy 


w^{z) = 


( 20 ) 


zm^{z) 7n^:(z) 

is analytic on C — , x+], verifies w^:{xf) = y/c^, and increases from cj^ to -foo when x increases 

from xt to -foo (see ||3l, section 3.1). Moreover, if (j)^{w) denotes function defined hy 

{w -f a‘^){w -f cr^c*) 


4>*{w) = 


w 


then, 0* increases from xf to -|-oo when w increases from a^^/cf to -foo. Finally, it holds that 

(l)^{w^{z)) = z 


( 21 ) 


( 22 ) 


for each z e C — [x^, , x+]. 

We denote hy Qn{z) and QAr(z) the so-called resolvent of matrices Z^rZ^ and Z^Z^v defined hy 
QAr(2^) = (ZatZ^ - , QAr(^) = (Z^Zat - z1nl)~^ 
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Then, the results of iH imply the following proposition. 

Proposition 1. • (i) The eigenvalue distribution of matrix converges almost surely towards the 

Marcenko-Pastur distribution or equivalently, for each z € C — IR^, 

-m^{z) ^Qa.s. (23) 

• (ii) For each e > 0, almost surely, for N large enough, all the eigenvalues of Z^rZ^ belong to 

[(T^(l — “ e, (7^(1 + + A ¥ c* < 1, and to — e, a^{l + + e] U {0} 

i/c* > 1. 

• (Hi) Moreover, j/otv, byr are (M—L+l)—dimensional deterministic vectors satisfying supjv(||aAr||, l|bAr||) < 
+00 , then it holds that for each z G C"*" 

ay (QAr(-s) - m*(z)I) by —)• 0 a.s. (24) 

Similarly, if ajvf and by are NL—dimensional deterministic vectors verifying supy(||ay ||, l|bAr||) < 

+ 00 , then for each z G C'*', it holds that 

ay (^Qy( 2 ) — m*(z)I^ by — )• 0 a.s. (25) 

Moreover, for each z G C+, it holds that 

ay (Qy(z)Zy) by —> 0 a.s. (26) 

Finally, for each e > 0, convergence properties d24l l25l \26\l hold uniformly w.r.t. z on each compact 
subset o/ C — [0, + e]. 

The proof is given in the Appendix. 

Remark 1. Proposition \I\ implies that in a certain sense, matrix ZyZy behaves as if the entries of 
7 jn were i.i.d because Proposition \I\ is known to hold for Ltd. matrices. In the i.i.d. case, m was 
established for the first time in /|9|/, the almost sure location of the eigenvalues of7iNX*p^ can be found 
in SB (see Theorem 5-11), while l l24l) . d25l) and < 1261) are trivial modifications of Lemma 5 of SB. 

We notice that the convergence towards the Marcenko-Pastur distrihution holds as soon as N ^ +oo 
and ^ c*. In particular, the convergence is still valid if iV = 0{M^) for each 0 < /3 < 1, or 

equivalently if L = 0 {M°^) for each 0 < a < 1. L can therefore converge towards oo much faster than 
N. However, the hypothesis that /3 > 1/3, which is also equivalent to L = 0{M°‘) with a < 2/3, is 
necessary to establish item (ii). 


March 30, 2015 


DRAFT 




11 


B. The K largest eigenvalues and eigenvectors o/XtvX^. 

While matrix Zj\f does not meet the conditions formulated in ||3l, Proposition [T] allows to use the 
approach used in ||3l, and to prove that the K largest eigenvalues and corresponding eigenvectors of 
XatX^. hehave as if the entries of Ztv were i.i.d. In particular, the following result holds. 


Theorem 1. 'We denote by s, 0 < s < K, the largest integer for which 

Xs ^ O' y/Tf 

Then, for k = 1,..., s, it holds that 

^ a.s. , ^ (Afc + 0-2)(Afc + 0-2c) , 

;; - > Pk — (PyXk) — -;- > x* . 


N^oo 


At 


(27) 


(28) 


Moreover, for k = s + ... ,K, it holds that 

Xk,N ^^^(1 + \/c7)^ (29) 

Finally, for all deterministic sequences of unit norm vectors (di tv). (d 2 ,Af). wo have for k = 1,..., s 

^l,N'^k,N'^k,N^‘2,N = 

/l*(/9fe)di TvUfc.AfUfc Tvd 2 ,Ar + o(l) a.s., (30) 

where function h^:{z) is defined by 


h,f{z) = 




( 31 ) 


w^{z){w^{z) + fj^c*) 

For the reader’s convenience, we provide in the appendix some insights on the approach developed in 
ll3i to prove (|2^ and (l29l) . For more details on (l30l) . see the proof of Theorem 2 in ||5l as well as the 
identity 


h^,{z) = 


zm^{zYm^{z) 

{zmYz)mYz))' 


where ' represents the derivation w.r.t. z. 


IV. Derivation of a consistent G-MUSIC method. 

We now use the results of section 17171 for matrix Xjv = /s/NL and Bjv = A(-^)(Sjv ® II). 

We recall that {Xk,N)k=i,...,M-L+i and (u^ Tv)fc=i,...,Af-i,+i represent the eigenvalues and eigenvectors 
of the empirical covariance matrix /NL, and that {Xk,N)k=i,...,K and (u^, Tv)fc=i,...,A' are the 

non zero eigenvalues and corresponding eigenvectors of (StvS^/AOIl) A(^)*. We recall that 

represents the orthogonal projection matrix onto the noise suhspace, i.e. the orthogonal complement 
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of the space generated by vectors {siM-L+i{Gk))k=i,...,K and that riN{ 6 ) is the corresponding MUSIC 
pseudo-spectrum 

m{0) = ^M-L+i{or HM-L+iio) 

Theorem [T] allows to generalize immediately the results of 151 and ifTSl concerning the consistency of 
G-MUSIC and MUSIC DoA estimators in the case L = 1. More precisely: 


Theorem 2. Assume that the K non zero eigenvalues (A^ 7 v)fc=i,...,it' converge towards deterministic 
terms Ai > A 2 > • • • > \k cind that 


Xk > (y\fc. 


Then, the estimator ^Ar(0) of the pseudo-spectrum r]]\f{6) defined by 

K 


VNiO) = i^M-L+liO)y 


E 


1 


=1 h lXk,N] 




(32) 


(33) 


verifies 


sup \fiN{0) - r]N{0)\ - ' 0, 

6»e[-7r,7r] 


(34) 


This result can be proved as Proposition 1 in |[5l. 

In order to obtain some insights on condition (l32l) and on the potential benefits of the spatial smoothing, 
we explicit the separation condition (l3^ when M and N converge towards 00 at the same rate, i.e. when 
^ d*, or equivalently that /3 = 1 and that L does not scale with N. In this case, it is clear that c* 

coincides with c* = d^,/L. It is easily seen that 


^a''-) ®II) Al'-l* ^{M-L+ 1/M) Am-l+1 ( 


^ • AlA, 


'-M-L-Vl 


(35) 


where • represents the Hadamard (i.e. element wise) product of matrices, and where B stands for the 

s s* 

complex conjugation operator of the elements of matrix B. If we assume that converges towards 

CJ C!* __ 

a diagonal matrix D when N increases, then 


• (AJ^Az,) converges towards the diagonal matrix 


D • Diag (A^Ai) = D. Therefore, 


N 


(AJ^Al) ~ D when is large enough. Using that 


M 


0 , 


we obtain that the separation condition is nearly equivalent to 


Xk (Am-l-i-iD a 


^M-L+l 


)> 


or to 




a'^y/df 

a'^y/df 


(36) 
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for each (M, N) large enough. If L = 1, the separation condition introduced in the context of (un¬ 
smoothed) G-MUSIC algorithms (Q) is of course recovered, i.e. 

Ax (A^AmD) > 


for each (M, N) large enough. If M is large and that L « M, matrix A^_^_|_^Am-l-i-i is close from 
A^Am and the separation condition is nearly equivalent to 


Ax (A^AmD) > 


O" 




Vl 


Therefore, it is seen that the use of the spatial smoothing scheme allows to reduce the threshold 
corresponding to G-MUSIC method without spatial smoothing hy the factor \/L. Hence, if M and 
N are the same order of magnitude, our asymptotic analysis allows to predict an improvement of the 
performance of the G-MUSIC methods based on spatial smoothing when L increases provided L « M. 
If L becomes too large, the above rough analysis is no more justified and the impact of the diminution 
of the number of antennas becomes dominant, and the performance tends to decrease. This analysis is 
sustained by the numerical simulations presented in section |Vl 


We define the DoA G-MUSIC SS estimates (0fc,x)fc=i,...,x by 

Ok,N = argmin \fjN{0) \ , (37) 

deik 

where Ik is a compact interval containing 6 k and such that n X/ = 0 for A; / 1. As in 151, (l34t as 
well as the particular structure of directional vectors aM-L-i-i(^) imply the following result which can 
be proved as Theorem 3 of 151 

Theorem 3. Under condition d32D . the DoA G-MUSIC SS estimates (0fc,x)fc=i,...,x verify 

M ~ 0 a.s. (38) 

for each k = 1,..., K. 

s s* 

Remark 2. We remark that under the extra assumption that converges towards a diagonal 

matrix, (13/ (see also m for more general matrices S) proved when L = 1 that (^k ,N - dk^ 

converges in distribution towards a Gaussian distribution. It would be interesting to generalize the results 
o/dH/ ([79|/ to the G-MUSIC estimators with spatial smoothing in the asymptotic regime AIR . This 
is a difficult task that is not within the scope of the present paper. 
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Theorem [T] also allows to generalize immediately the results of lITSl concerning the consistency of the 
traditional estimates in the case L = 1. In particular, while the traditional estimate of 

the pseudo-spectrum is not consistent, it is shown in lITSl that if L = 1, then the arguments of its local 
minima are consistent and verify 

M -Gk)^0 a.s. (39) 

for each A: = 1,... , if if the separation condition is verified. The reader can check that Theorem [U allows 
to generalize immediately this behaviour to the traditional DoA MUSIC estimates with spatial smoothing 
in regime More precisely, the following result holds. 

Theorem 4. Under condition f U2l) . the DoA traditional MUSIC SS estimates verify 

M -Gk)^0 a.s. (40) 

for each k = 1,..., K. 

Remark 3. It is established in / li^l/ in the case L = 1 that if converges towards a diagonal matrix, 
then {^k \ ~ ^ Gaussian behaviour, and that the corresponding variance coincides with 

the asymptotic variance of {Ok,N — In particular, if L = 1, the asymptotic performance of 
MUSIC and G-MUSIC estimators coincide. It would be interesting to check whether this result still holds 
true for the MUSIC and G-MUSIC estimators with spatial smoothing. 

Theorems |2] and |3] as well as (l39l ) assume that the DoAs {9k)k=i,...,K are fixed paramefers, i.e. do nof 
scale wifh M. Therefore, fhe rafio 

miiifc^; \9k-0i\ 

{2-k)IM 

converges fowards +oo. In practice, fhis confexf is fhus able fo model pracfical sifuafions in which 
sup;.^; \ 9k — 9i\ is significanfly larger fhan fhe aperture of fhe array. In fhe case L = 1, ifTSl also 
addressed fhe case where fhe DoA’s {9k^N)k=i,...,K depend on N,M and verify ~ 

This confexf allows fo capfure pracfical sifuafions in which fhe DoA’s are spaced of fhe order of a 
beamwidfh. In order fo simplify fhe calculations, ifTSl considered fhe case K = 2, 02 ,n = ^i, 7 V + ^ and 
where mafrix —)• I 2 . However, fhe resulfs can be generalized easily fo more general sifuafions. If 

is shown in ifTSl fhaf fhe G-MUSIC esfimafes still verifiy (l40l ). buf fhaf, in general, M — 9]^ does 
nof converge fowards 0. The resulfs of ifTSl can be generalized immediafely fo fhe confexf of G-MUSIC 
estimators wifh spatial smoofhing in regime (fTTI) . For fhis, we have fo assume fhaf 02,7V = ^i, 7 V + 
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(in ifTSl . M and N are of the same order of magnitude so that the assumptions 02,7V = Oi,N + ^ and 
02,7V = 01,7V + are equivalent), and to follow the arguments of section 4 in ifTSl . The conclusion of 
this discussion is the following Theorem. 

s s* 

Theorem 5. Assume K = 2, 02,at = 0i,7v + and that —)• I 2 . If the separation condition 

1 — |sincK/2| > (T^c* (41) 

holds, then the G-MUSIC SS estimates (0fc,7v)fc=i,2 defined by 

Ok,N = argmin |j77v(6')| , (42) 

where Ik,N = [^fc,7V ~ Sk,N + far e small enough, verify 

M (^Ok^N ~ ^k,N^ 0 a.s. (43) 

In general, the traditional MUSIC SS estimates defined by d42l) when the G-MUSIC estimate i}7v(0) ^ 
replaced by the traditional spectrum estimate {0) are such that M — 0fc,7V^ does not converge 
towards 0. 


V. Numerical examples 

In this section, we provide numerical simulations illustrating the results given in the previous sections. 
We first consider 2 closely spaced sources with DoAs 0i = 0 and 02 = and we assume that 
M = 160 and N = 20. The 2 x N signal matrix is obtained hy normalizing a realization of a random 
matrix with A/’c(0,1) i.i.d. entries in such a way that the 2 source signals have power 1. The signal to 
noise ratio is thus equal to SNR = Xja^. Table U provides the minimum value of SNR for which the 
separation condition, in its finite length version (i.e. when the limits {Xk)k=i,...,K and c* are replaced by 
{Xk,N)k=i,...,K and Cat respectively) holds, i.e. 

( tt ^)-! = ^^{M-L + 1)/NL 
^K,N 

It is seen that the minimal SNR first decreases but that it increases if L is large enough. This confirms 
the discussion of the previous section on the effect of L on the separation condition. 

In figure [T] we represent the mean-square errors of the G-MUSIC SS estimator 0i for L = 2,4, 8,16 
versus SNR. The corresponding Cramer-Rao bounds is also represented. As expected, it is seen that the 
performance tends to increase with L until L = 16. In figure |2j L is equal fo 16, 32, 64, 96, 128. 

For L = 32, it is seen that the MSB tends to degrade at high SNR w.r.t. L = 16, while the performance 
severely degrades for larger values of L. 
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L 

2 

4 

8 

16 

32 

64 

96 

128 

SNR 

33.46 

30.30 

27.46 

25.31 

24.70 

28.25 

36.11 

51.52 


TABLE I: Minimum value of SNR for separation condition 



Fig. 1: Empirical MSE of G-MUSIC SS estimator 6 i versus SNR 



Fig. 2: Empirical MSE of G-MUSIC SS estimator 9i versus SNR 
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In Figure [3l parameter L is equal to 16. We compare the performance of G-MUSIC SS with the 
standard MUSIC method with spatial smoothing. We also represent the MSE provided hy G-MUSIC and 
MUSIC for L = 1. The standard unsmoothed MUSIC method of course completely fails, while the use 
of the G-MUSIC SS provides a clear improvement of the performance w.r.t. MUSIC SS and unsmoothed 
G-MUSIC. 



Fig. 3: Empirical MSE of different estimators of 9i when E=16 


We finally consider the case L = 128, and compare as above G-MUSIC SS, MUSIC SS, unsmoothed 
G-MUSIC and unsmoothed MUSIC. G-MUSIC SS completely fails because L and M are of the same 
order of magnitude. Theorem |2] is thus no more valid, and the pseudo-spectrum estimate is not consistent. 

We now consider 2 widely spaced sources with DoAs 6 i = 0 and 62 = 5^, and keep the same 
parameters as above. We consider the case L = 16, and represent in Fig. [5] the performance of MUSIC, 
G-MUSIC, MUSIC-SS, and G-MUSIC-SS. It is first observed that, in contrast with the case of closely 
spaced DoAs, MUSIC-SS and G-MUSIC-SS have the same performance when the SNR is above the 
threshold 6 dB. This is in accordance with Theorem |4l and tends to indicate that, as in the case L = 1, 
if converges towards a diagonal matrix, then the asymptotic performance of G-MUSIC-SS and 

MUSIC-SS coincide (see Remark 1^. The comparison between the methods with and without spatial 
smoothing also confirm that the use of spatial smoothing schemes allow to improve the performance. 
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Fig. 4: Empirical MSE of different estimators of 6 i when E=128 



Eig. 5: Empirical MSE of different estimators of 0i when E=16 and widely spaced DoAs 


VI. Conclusion 

In this paper, we have addressed the behaviour of suhspace DoA estimators based on spatial smoothing 
in asymptotic regimes where M and NL converge towards oo at the same rate. Eor this, we have evaluated 
the behaviour of the largest singular values and corresponding singular vectors of large random matrices 
defined as additive low rank perturbations of certain random block-Hankel matrices, and established that 
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they behave as if the entries of the block-Hankel matrices were i.i.d. Starting from this result, we have 
shown that it is possible to generalize the G-estimators introduced in IH, and have deduced from lITSl 
their properties. 


Appendix A 

Insights on the proof of (l2^ and (|2^ . 


We first recall that Q established that for 1 < A; < AT, if Xk,N does not converge towards a limit 
strictly greater than xf, then Xk,N converges towards xf. We have therefore to evaluate the behaviour 
of the eigenvalues of that are greater than + e for some e > 0. 

If C represents a P x Q matrix, we denote by C the {P + Q) x {P + Q) hermitian matrix defined by 


C = 


0 c 

C* 0 


(44) 


Then, fhe non zero eigenvalues of C coincide wifh fhe (posifive and negafive) square roofs of fhe 
non zero eigenvalues of mafrix CC*, and fhe corresponding eigenvectors are fhe (P + (5)-dimensional 
vectors (a^, ±b^)^ where (a^, b^) represenf fhe pairs of leff and righf singular vectors of C. Therefore, 
A > 3:+ + e is eigenvalue of XatX^ if and only if y/X > (x+ + is eigenvalue of matrix X^. We 
consider the singular value decomposition 


Bn = 


of matrix B n and express X^y as 


which can be written as 

0 

Zjv* 


0 

Zn 

+ 

0 


v* 

0 

0 







Ztv 

+ 

Bn 

0 

J 

1 

o 

_1 

0 

0 





D D* 


where J is defined by 
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Consider x > {x^ + Then, x is not a singular value of ZatZ^, and therefore, not an eigenvalue 
of Z^r. Therefore, it holds that 

det(XjY — a:I) = det(Z^ — xl + DJD*) 

= det(Z^-xI)det(I + (Zjv-xI)~^DJD*) 

= det(Z^ - xI)det(l 2 K + JD*(Z^ - xI)"^D) 

after noticing J = J”^. For to G C — [—(x+ + (x+ + we denote hy SAr(to) the 2K x 2K 

matrix defined hy 

StvM = l2i^ + JD*(Z^ - tt;I)-^D 


Using the identity 

(Z^ - tol)-' 


we obtain immediately that 


toQAr(to^) Qn{w'^)’Z‘N 
ZlfQNiw"^) wQn{w‘^) 


= 1k + A]i^\J%Z%QN{w^)VN 

(S7v(w;))i,2 = wA]i^\J*j^QMiw^)'(JNA]i^ 

(Sjv(^/^))2,1 = u;U:^Qiv(w;")U,v 

(S^(tt;))2 2 = Ik + lJ*NQN{w^)ZN'(JNA]i^ 

Item (iii) of Proposition [T] implies that the elements of Sat (to) converge almost surely, uniformly on the 
compact subsets of C — [—(x;)" + {x^ + e)^/^] towards the elements of matrix S*(to) defined by 

Ia" wfh^{w‘^ )A 
wm{wP‘)lK liv 

If is easy fo check fhaf det(S 7 v('»^)) and det(S*(to)) are funcfions of to^. We define funcfions sn and s* 
on C — [0,x+ + e] by sn{w^) = det(SAr(m)) and s^{w‘^) = det(S*(to)). If is clear fhaf almosf surely, 
sn{z) s^{z) uniformly on fhe compacf subsefs of C — [0, (x+ + e)]. Therefore, in order fo precise fhe 
behaviour of fhe eigenvalues of XatX^ fhaf are greater fhan x+ + e (i.e. fhe solufions of fhe equation 
sj\f{x) = 0 greater fhan x^ + e), if is firsl useful fo characterize fhe solutions of fhe equation s*(x) = 0. 
The equafion s* (x) = 0 is equivalenf fo 

^k=i (1 “ Afcxm*(x)?7i*(x)) = 0 


S*(to) = 
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or equivalently to 


w^{x) = Xk 


for k = 1,... ,K. Using the properties of function m*, we obtain immediately that if e < ps — xf = 
0*(As) — x^, then the solutions of s*(x) = 0 that are greater than xf + e coincide with the {pk)k=i,...,s 
defined by pk = (j)*{Xk) for k = 1,..., s. Using this, it can be proved using appropriate arguments that, 

ataos, surely, for large enough, .hen .he s greats, eigenvalues (Aa,„)a..., of are grea.er 

than xf + e, and that Xk,N —s- Pfc for A: = 1,..., s [J. 


Appendix B 

Proof of Proposition [T] 

The proof of Proposition [T] is based on the results of lEl. In order to explain this, we denote by W^r 
the NL X [M — L + 1) matrix defined by 

The variance of the entries of Wat is equal to Therefore, matrix Wat is similar to the matrices 

studied in HI except that the integers (M, N) in HI should be exchanged by {N, M — L + 1). In particular, 
after this replacement, it is clear that the asymptotic regime (fTTl) coincides with the regime in HI- In order 
to recall the results of HI. we denote by tN{z), iN{z), t^{z) and i^{z) the Stieltjes transforms of the 
Marcenko-Pastur distributions of parameters (cj^,c^^), cat), and (cr^cT^,c*). Moreover, 

Qn,w{z) and Qn,w{^) represent the resolvents of matrices WatW^^t and WJ^Wat respectively. It is 
shown in HI (see Section 6) that the eigenvalue distribution of WatW^^t converges almost surely towards 
c-D a statement equivalent to 

-^Tr(QAr,iv( 2 ;)) - Uiz) 0 a.s. 

or to 

- U{z) ^ 0 a.s. 

for each z e C+. As we have 

Z%Zn = cnWnW*^ ( 45 ) 

'The arguments used in m require the uniform convergence of sjv towards s* on the set Re(z) > xif + e, a property that is 
not established in Proposition [T] However, the proof of the contuinity lemma 2.1 in m can be simplihed, and only needs the 
uniform convergence on compact sets. 
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it holds that the resolvent Qn{z) of is equal to 

QAf(^) = C^^Qn,w{zc'^^) 

As cn —^ c*, it is clear j^Tr(Q 7 v(- 2 )) behaves as Similarly, ]g:rc+iTr(QAr( 2 :)) 

behaves as ^ jj_\^_)_^ Tr(Qjv,w(-gc7^)). Therefore, for each z G C+, it holds that 

■^Tv{Qn{z)) - c~^ U{zc~^) 0 a.s. 

and that 

^^Tr(QAr( 2 ;)) - c~^ i4zc~^) 0 a.s. 

Using Equations (fT^fT^ . it is easy to verify that m^{z) = c~^ i^:(zc~^) and m^{z) = c~^ t^{zc~^). This 
establishes (|2^ and the convergence of the eigenvalue distribution of ZatZ^ towards 

Asymptotic regime (fTTl) implies that L = 0{M°^) = 0{{M — L + 1)“) where a < 2/3. Therefore, ||8l 
implies that for each e > 0, almost surely, for N large enough, the eigenvalues of WatWI^ are located 
in [(T^(l — \/+ y/cT^)^] U {0}l(cT^ > 1). (1^ and the convergence of cn towards c* 
lead immediately to item (ii) of Proposition [T] 

Using the same arguments as above, (l25l) appears as a consequence of 

(QAT.tyC-^) - U('^)I) ijAf0 o.s. (46) 

While (|4^ does not appear explicitely in IH, it can be deduced rather easily from the various intermediate 
results proved in ||8l. For this, we first remark that 

^*N {^N,w{z) - t^{z)T) Ijat = (QAr,w( 2 ;) “ E(QAr,vv(^))) byr + (E(QAr,vy( 2 ^)) - U{z)l) byr 

and establish that the 2 terms at the right hand side of the above equation converge towards 0. In order to 
simplify the notations, we denote by ^ the first term. The almost sure convergence of ^ towards 0 follows 
from the Poincare-Nash inequality (see e.g. Proposition 2 of ||8l). Exchanging (M, N) by {N, M —L + 1) 
in Proposition 6 of HI, we obtain immediately that E|^p = 0( = 0{^). As L/M —)• 0, this 

implies that ^ converges in probability towards 0. In order to prove the almost sure convergence, we 
briefly justify that for each n, it holds that 

E|^|2^ = C)((L/M)”) (47) 

Wh can be established by induction on n. As mentioned above, WT\ is verified for n = 1. We now 
assume that it holds until integer n — 1, and prove (1471 ). For this, we use the obvious relation: 

E|C|2" = (E|er)VVar(n 
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Using the Poincare-Nash inequality as in the proof of Proposition 6 of HI, we obtain easily that 

Var(n<C'^E(|e|2-2) 

where C is a constant that depends on 2 ; hut not on the dimensions L, M, N. As Wt\ is assumed to hold 
until integer n — 1, this implies that Var(,^”) = O ((L/M)"^). The Schwartz inequality leads immediately 
to 

(Elerf <E(|enE(|Cp”-2) 


which is a O {{L/M)'^) term. This establishes (l47l) . As L = 
(L/M)^ verifies 


(L/M)^ 


1 

j\,fl+2-3a 


0{M^) with a < 2/3, it is clear that 


Therefore, WT\ for n = 3 leads to 

As 2 — 3 q; > 0, the use of the Markov inequality and of the Borel-Cantelli lemma imply that ^ converges 
towards 0 almost surely as expected. 

It remains to justify that (E(Q 7 v,iy ( 2 )) — f*( 2 ;)I) hj\f converges towards 0. Although it is not stated 
explicitely in m, it can immediately deduced from Eq. (5.3) in Proposition 8, as well as on Corollary 
1, Theorem 2, and formula (7.3). 


(l24l) is equivalent to 


[QN,wiz) - U{z)l) Uat 0 a.s. 


It can be proved as above that 


( QN,wiz) - E(QAr,w(^))) bAT —0 a.s 


(48) 


and establish that 

(e(Q,v,iv(^)) - U(^)l) ^ 0 (49) 

for each 2 ; G C+. The behaviour of matrix E(QAr^vv(-2)) is not studied in iHl. However, it can be evaluated 
using the results of m. For this, we first simplify the notations and denote by W, W, Q, Q the matrices 
Wtv, Wat, Q 7 v,ty( 2 :), and QN,wiz). Moreover, Q is a NL x NL block matrix, so that we denote by 
its entry (ii + (m - l)L,i 2 + (n 2 - 1)L). 
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As in iH, we denote by and ^i-) the operators defined by 


'M-L+1,L\ 
-{N) 


1 




NL 


L-l 




-W /r,\ - T-Wj 

i=-iL-l) 

where Jl is the L x L ’’shift” matrix defined by (J L)ij = 6{j — i = 1). 

Replacing integers (M, N) by integers {N, M — L + 1) in Equation (4.6) of |0, we obtain that 


E 




(7 \ 

- / , ll<fc-i<M-L+lE 

i=-(L-l) 


ii,k^ 
2 L-l 


J,*2 


a 


M- L + 




r('^)(Q)(*)(QW)"E_,(W^ 


*\n2 


with l<j,k<M — L + 1, 1 < ni,n 2 < N, I < ii,i 2 < L- 

Setting u = k — i, the second term of the righthandside of the above equation can also be written as 

2 M —Z/+1 

— J] E[rW(Q)(fc-u)l_(i_i)<fc_„<z._i(QW)-^jW*)-;„; 


U=1 


Now setting n = ni = n 2 , f = ii = ^ 2 ^ and summing over n and i, we obtain 

2 

E(W*QW)._, = 


a 


CN 


^M-L+l 


E J] TW(Q)(A:-n)l_(z._i)<fc_„<i_i(W*QW),-, 


U=1 


and using that tI"I(Q)(J; - we get that 

E(WQW),,, = £ (ri,"Ai,L(E(Q))) 

o 

We express matrix Q = E(Q) + Q, and obtain that 


a 


—E r 


k,j CN 


-(TV) 


(Q)W^Q^W 


' M-L+1,L 


k,j 


E(WQW)_,., = £ (rrA,,i(E(Q)))^ , 


o 

CN 


rifJ^+i,L(E(Q))E(W^Q^W)) 
2 


a 


' k,j 


—E r 


CN 


(N) 


(Q)W^Q^W ) 


' M-L+l,L 


' k,j 


Noticing the equation, 


W^Q^W = W 


T\■\TT^ 


we obtain that 


E (WQW) = (E(Qq)-^E (wWq) (E(q^))-^E WWQri'!.Ei,i(Q) 


CN 


CN 


CN 
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Moreover we notice that 


Therefore, it holds that 


W*QW = QW*W = W*WQ = I + zQ 


I + zE(Q) = JE(Q^)) - —(I + (E(Q^)) + A 

Cjv ctv 


where 


a 


A = -^E ( WWQ4''_T,,i(Q) 

C_/V 


This leads to the equation 


zE{Q) I + —)) = _i + A 

CN 


Lemma 2 of 181 (used when (M, N) is replaced by (M — L + 1, N)) implies that matrix 

CN 

is invertible for z e C'*', and that its inverse, denoted H, verifies 


(50) 


(51) 


(52) 


IHII < 


Im(z) 


for 2 ; G C"''. (l52l) implies that 

Therefore, (|4^ is equivalent to 


TTi 

E(Q) =-+ AH 


T 


H 


T 


(53) 


—— - Uiz)l + AH" J biv ^ 0 

Using the same technics as in Proposition 8 (see Eq. 5.3) of || 8 l as well as (l5^ . we obtain immediately 
that 


It thus remains to establish that 


a;^AH" bjv ^ 0 


a^ (-] bjv —y 0 


(54) 


For this, we use the identity 


H" 


-U(z)I = -H'^ 


+ (H^)-Mt;(z) 


zt^{z) 


t^:{z) and U(z) satisfy the relation 7 ^ 77 ^ = 1 + Hence, the right hand side of the above equation 


can be written as 


a 


zi^(z) 
2 


-H" (-1-U(z))I + I 

CN 


^<-L+i^L (E(Q^))) i.(^) 
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Corollary 1 and Theorem 2 of HI imply that 

{E(Q^) - 

if z € C+. This and dS^ leads to 

and to (l54l) . This completes the proof of (l49l) . 


We now establish (l26l ). For this, we first remark that for each 0 G IR, the distribution of matrix Z^re*^ 
coincides with the distribution of Ztv- Therefore, it holds that 

E (Q,v(^)Zjve''') =E{Qn{z)Zn) 

which implies that E {QN{z)Z]\f) = 0. In order to complete the proof of (|2^ . it is sufficient to establish 
that if we denote by kn the random variable kn = (Q 7 v( 2 ;)Z 7 v) hj\f, then, for each p > 1, it holds 

that 

= (55) 

Choosing p large enough leads to kn — E(k; 7 v) = kn —^ 0 a.s. as expected. (1551) can be proved as above 
by using the Poincare-Nash inequality. 

We finally justify that for each e > 0, (l24l l25l f26\i hold uniformly w.r.t. z on each compact subset of 
C — [0, xf + e]. We just prove that it the case for (|2^ . By item (ii), almost surely, function z ^ kn{z) 
is analytic on C — [0,a;+ + e]. We use a standard argument based on Montel’s theorem ( lfT2]| . p.282). We 
first justify that for each compact subset /C C C — [0, + e], then it exists a constant r] such that 

sup \ KNiz)\ < T] (56) 

z&K 

for each N large enough. We consider the singular value decomposition of matrix Zat : 


ZiV = TatAat©^ 

where Atv represents the diagonal matrix of non zero singular values of Zjv- kn{z) can be written as 

kn{z) = (A^ - ziy^ An&*n^n 

Therefore, it holds that 


\>iN{z)\ < 




-Zl) "An 


laArll lib 


'N\ 
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Item (ii) implies that the entries of are located into [0, + e]. for each N large enough. Therefore, 

for each 2 ; G /C, it holds that 




-. 1 )' 




< 


dist([0, x* + e],/C) 

The conclusion follows from the hypothesis that vectors a^v and b^r satisfy sup^(||aAr||, ||bAr||) < + 00 . 
(l56l) implies that the sequence of analytic functions (KAr) 7 v>i is a normal family . Therefore, it exists a 
subsequence extracted from {kn)n>i that converges uniformly on each compact subset of C — [0, + e] 

towards a certain analytic function k*. As (l26l ) holds for each 2 G C+, function k* is identically zero. We 
have thus shown that each converging subsequence extracted from {kn)n>i converges uniformly towards 
0 on each compact subset of C — [0,3:+ + e]. This, in turn, shows that the whole sequence converges 
uniformly on each compact subset of C — [0, x+ + e] as expected. 
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